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We measure the penetration depth Xab{T) in Ba(Fei_a;Coa;)2As2 using local techniques that do 
not average over the sample. The superfluid density psiT) = l/Xab{T)^ has three main features. 
First, Ps(T ~ 0) falls sharply on the underdoped side of the dome. Second, \ab{T) is flat at low T at 
optimal doping, indicating fully gapped superconductivity, but varies more strongly in underdoped 
and overdoped samples, consistent with either a power law or a small second gap. Third, psiT) 
varies steeply near Tc for optimal and underdoping. These observations are consistent with an 
interplay between magnetic and superconducting phases. 
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The magnetic penetration depth A is one of the most 
important length scales in superconductors [I]. Its tem- 
perature evolution is a sensitive probe of the supercon- 
ducting gap structure [21 [3] and its value is related to 
the density of electrons in the superconducting state, 
Ps{T) = 1/A(r)^. Comparing Ps{T) from samples across 
the phase diagram of a family of iron pnictide supercon- 
ductors can shed light on the pairing mechanism as 
well as on the relation between superconductivity and 
adjacent phases [5]. The family of Ba(Fei_2;Coa;)2As2 is 
a good candidate for such studies because single crys- 
tals can be grown cleanly with relatively sharp super- 
conducting transitions [BJ [7] , because the magnetic and 
structural transitions are well characterized [5HTU]. and 
because other careful studies of quasi-particle excitation 
measurements such as thermal conductivity [11] [12] and 
Raman scattering |13| have been successfully performed 
across the phase diagram. 

A is difficult to measure accurately. In the pnictides, 
the proximity of the magnetic phase to superconductiv- 
ity on the underdoped (UD) side of the superconduct- 
ing dome prevents measurement of A by ^SR [Hj. Bulk 
measurements by microwave- and RF-based techniques 
are made difficult by complex sample topography and 
inhomogeneity, which can explain the significantly dif- 
ferent results among nominally similar samples of both 
AA = A(r) - A(0) [BHU] and A [IHj measurements. 
Local-probe studies of ps (T) across the doping range are 
strongly desirable, because such measurements can ob- 
tain A even when the magnetic order is adjacent or co- 
existing [21) . and can reduce the error from topography 
and inhomogeneity [22 • 

In this Letter, we report local measurements of XabiT), 
the penetration depth for screening currents flowing in 
the ab plane in a set of Ba(Fei_a;Coa;)2As2 single crys- 
tals grown from self-flux [7|. We measure at Co compo- 
sition X across the superconducting dome: underdoped 
(UD) X = 0.045, 0.049, 0.051; optimally doped (OptD) 



X = 0.07; and overdoped (OD) x = 0.085, 0.11 ^. Our 
measurements average over a few microns. The positions 
we choose to measure show strong, uniform diamagnetic 
response, and are at least 15 microns away from topo- 
graphic steps larger than 0.5 micron. We are able to re- 
solve well-formed vortices to rule out granularity on sub- 
micron scales except at x = 4.5%. We observe system- 
atic evolution of Ps{T) with x. Our observations suggest 
strong correlation between magnetism and superconduc- 
tivity. 

We use magnetic force microscopy (MFM) to measure 
\ab{T) and AAa6(T) = \ab{T) - Aab(O) irrespective of 
adjacent magnetic states from 5 K to the superconduct- 
ing transition temperature Tc pi]. We also use scanning 
Superconducting QUantum Interference Device suscep- 
tometry (SSS) in a '^He cryostat to measure AAaf,(r) 
from 0.4 K to 7 K [22]. Both techniques measure the 
diamagnetic response to a local field source (either the 
MFM tip or the SQUID field coil) in the Meissner state. 
This is measured by SSS through the mutual inductance 
between the field coil and the pickup loop, and by MFM 
through the derivative of the vertical component of the 
force between the magnetic tip and the sample dF^/dz, 
where z is the tip-sample distance. The diamagnetic re- 
sponse can be approximated as a time-reversed mirror of 
the source reflected about a plane Xab below the surface 
when the field source is much further than Xab above the 
surface 23J. Under this approximation, changes in z and 
A are equivalent, allowing model independent measure- 
ment of AAab(r). By properly modeling the MFM tip- 
superconductor interaction, we also obtain Xab{T). More 
details on the techniques are provided in previous publi- 
cations [2X1122]. 

We determine Xab{T) from 5 K to Tc by MFM within 
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15% error, mostly from uncertainty in the MFM tip ge- 
ometry. We use the same MFM tip at a; = 4.5%, 7% 
and 11%, reducing the relative error of Xab{T) among 
these samples. AAa&(T) by both techniques has 7% 
error, mostly from calibration uncertainty of the scan- 
ner. At most X, we measure at least two samples from 
the same growth batch, one by MFM and one by SSS. 
When AXab{T) measured by both the techniques overlaps 
in the common temperature range, we offset AAQf,(T) 
from SSS by XabiT = 5K) from MFM, and obtain 
Ps{T) = l/\ab{TY over the full temperature range, as 
shown in Fig.[l] However, /S.\ab{T) measured by the two 
techniques do not agree at a; = 4.9%, presumably due 
to sample variations. At a; = 11%, we observe different 
AAab(T) by SSS at three locations separated by hundreds 
of microns (shown in Fig. |3]) , one of which matches the 
MFM measurement. 

We observe a systematic change of ps (T) with Co dop- 
ing X across the superconducting dome that can be char- 
acterized by three trends. First, the zero temperature 
value Ps(0) (Fig. [2]) peaks at OptD. It is strongly reduced 
in the UD regime, falling sharply as the magnetic order 
onsets. For UD samples, ps{0) falls more quickly than 
Tc, while on the OD side, Ps{0) falls less rapidly than Tc. 
This observation is different from a previous measure- 
ment by bulk technique that reported PsiO) increasing 
with doping across the doping dome [5D]. 

Second, the magnitude of AAa6(T) [Fig. [ij:] at low T 
increases away from OptD on both sides of the dome. At 
X — 7%, AAab(T) remains flat within 0.5 nm below 3 
K, indicating fully gapped superconductivity, consistent 
with the proposed order parameter s± [31. We use a 
phenomenological single-gap fit to Ps{T) at a; = 7%, with 
the gap function [3] 

A(T) = A(0)tanh(||yr(|r;)) (1) 

where a is a free parameter that characterizes the rise 
of ps(T) below Tc. This fuh single gap fit gives A(0) = 
2.0/cbTc and a = 1.7, and, as shown in figure la, ade- 
quately describes the measured temperature dependence 
for a; = 7%. 

Due to the steeper AAab(r) away from OptD, the 
full single gap fit does not work at other dopings. The 
low temperature variation can be described by a power 
law, but the full temperature dependence is also well de- 
scribed by a two-gap fit with one gap being small. In 
cases where the SSS and MFM data do not agree or 
where only one is available, only one of the fits is pos- 
sible. We use a phenomenological two-gap fit, Ps{T) — 
ppi{T) + (1 - p)p2{T) (Eq. 2), where pi, p2 are given 
by Eq. 1 with gaps Ai, A2. We fix Tc of the two gaps 
to be identical and set a for the smaller gap to be 1 to 
minimize the number of free parameters. Best fit val- 
ues are shown in Table [l] In the power law description. 
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FIG. 1; The temperature dependence of ps{T) = l/\ab{Tf 
and AAai,(T) shows three systematic trends, a: PsiT) vs 
T. Solid line: a single-gap equation (Eq. 1) fits only x — 7%. 
Dashed lines: phenomenological two-gap fits (Eq. 2) work well 
at every doping. Fit parameters in Table. At a: = 8.5%, we 
did not measure Xab{T), so we offset A\ab[T) by the mean of 
Xab{0) at a; = 7% and x = 11%. b: ps{T)/ps{0) vs T/Tc. At 
X = 5.1%(UD) and x — 7%(OptD), ps(T) rises more sharply 
than for MgB2 from Ref. [24], single-band weakly coupled 
BCS theory, or a; = ll%(OD). c: A\ab{T) down to 0.4 K 
measured by SSS at the indicated x. Dashed lines: power-law 
fits (Eq. 3) describe the data well up to O.STc with the power 
n fixed at 2.5 (fit parameters in Table. The amplitude 
A increases away from optimal doping, where AAai,(r) is so 
fiat as to be consistent with exponential behavior. Successive 
data sets are offset vertically by 5 nm for clarity. /S.\ab(T) at 
X = 4.9% measured on a different sample by MFM at T > 5 K 
is also plotted. 



AAafc(r) ^T" (Eq. 3). The best fit value of n varied 
from 2.3 to 2.9, but the error bars for all data sets are 
consistent with n = 2.5. We therefore fixed n = 2.5 for 
simplicity. As shown in Table |l) the fitted magnitude of 
A increases away from OptD. 

The third trend we observe is that near Tc, Ps{T) of the 
OptD and UD rises faster with decreasing temperature 
than weak-coupling BCS theory, MgB2 or OD (Fig. [T]3). 
This feature is characterized in the full-gap fits: a at a; = 



3 



E 400 



,300 

200 

60 

40 
20 




10 ^ 




A 





□ Tc(bulk) 






♦ Tc (local) 












VTsm 




n 



4 6 8 10 
Co doping x (%) 



12 



FIG. 2: a: Aab(O) and pa(0) as extrapolated from Fig[lja), 
showing the rapid drop of Ps(0) on the underdoped side. For 
X = 4.5%, pa{5K) is plotted because the extrapolation is un- 
clear, b: Temperatures of structural, magnetic and supercon- 
ducting transitions Tstru, Tafm and Tc(bulk), from resistiv- 
ity measurements of samples from the same growth as studied 
here, rc(local) is from linear extrapolation of pa{T) to zero 
from the MFM data presented here. 



TABLE I: Fit parameters for curves in Fig.[l]and|3} for the 
two-gap (Eq. 2) and power-law (Eq. 3) models. Power law fits 
are based on SSS data and two-gap fits on MFM or combined 
MFM and SSS data. Where both power-law and two-gap fits 
are shown, SSS and MFM data on different samples agreed. 
M. X = 4.9%, 5.1% and 7%, AAai,(r) at muhiple positions 
separated by at least 100 microns are identical. For x=ll%, 
SSS results from 3 positions are shown. Tc measured on the 
same sample by MFM or SSS are also provided. 
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(nm/K^-^) 



4.9%-l 

4.9%-2 

5.1% 

7% 

8.5% 

11% 

11% 

11% 



15.8 
15.5 
18.6 
22.4 
19.6 
13.5 



2.2 0.6 



2.3 
3.3 
1.9 
1.7 



0.7 
1.3 
0.6 
0.6 



1.6 0.90 

0.26 (3 positions) 

1.6 0.90 0.09 (4 positions) 

1.7 0.70 0.02 (4 positions) 

1 0.92 

0.9 0.87 0.38 (position 1) 

0.27 (position 2) 
0.18 (position 3) 



5.1% and 7% are significantly larger than the BCS value 
a = 1, suggesting a steeper rise oi ps{T) with cooling than 
the weakly-coupled BCS would give. The same feature 
is shown in ps(T)/ ps{0) vs T/Tc plots when comparing 
with other superconductors. The curves at x = 5.1% 
(UD) and 7% (OptD) have similar slopes near T^,. Both 
are steeper than that of MgB2, a weakly coupled two-gap 
BCS superconductor; or a; = 11% (OD). 

Local probes allow us to measure at multiple positions 
on multiple samples and to examine real-space inhomo- 
geneity. On all samples we measured with MFM, we 



obtain identical Xab{T) and AXab{T) within errors on 
positions several microns apart. On all samples studied 
with SSS, we measured three or four positions separated 
by hundreds of microns. For the x = 4.9%, 5.1%, and 
7% samples, the AAob(T) measurements were identical 
within errors at each location. We observe spatial vari- 
ation in AXab(T) on one x — 11% sample as shown in 
Fig. [3] resulting in the variations in fit parameters shown 
in Table|Tj In principle, such spatial variations could come 
from doping inhomogeneity in the sample, or from rough- 
ness of the sample surface, or from other variations. In 
order to quantify doping inhomogeneity in the sample, 
we performed x-ray microanalysis after the SSS measure- 
ment was complete. The 3 fim beam was scanned across 
a 1.2 mm long line in 20 fj,m steps. We observe no sys- 
tematic change in the doping level to within ±2% of the 
total Co concentration. We do not think the spatial vari- 
ation comes from sample surface roughness because we 
checked the topography by susceptometry scans in-situ 
and by optical microscopy after the measurement, and 
because roughness is expected to overestimate AA, while 
the three values are either smaller than or the same as 
the MFM result. 

The systematic trends that we have observed in Ps{T) 
should be considered in light of properties in the super- 
conducting states that are expected to evolve with dop- 
ing, including the structure of the gap in fc-space on mul- 
tiple bands [13l [25] , magnetic scattering, other forms of 
scattering [26], and transfer of spectral weight to spin 
fluctuations and the magnetic phase [57]. In particular, 
the three trends of Ps{T) can be accounted for by the 
interplay between magnetism and superconductivity. 

The first trend, ps (0) dropping more rapidly on moving 
towards UD than towards OD, follows naturally from the 
fact that the structural and magnetic transitions lead to 
significant Fermi surface reconstruction [13 HH], resulting 
in smaller electron and hole pockets and therefore fewer 
charge carriers for the superconducting state. In the UD 
cuprates, the reduction of ps{Q) approximately following 
Tc 29] has been often attributed to phase fluctuations 
of the superconducting state [3D]. We observed an even 
faster drop of ps(0) than Tc, consistent with the scenario 
that coexisting order, e.g. magnetic order, removes a 
large number of itinerant carriers that might otherwise 
enter the superconducting condensate. 

The second trend, weakening of fully-gapped behavior 
away from OptD, agrees with heat transport measure- 
ments which have also reported an increase in low-energy 
quasi-particle excitation on either side of OptD [H] . The 
observation indicates strong pair-breaking scattering or 
anisotropic superconducting gap structure in the s± pair- 
ing symmetry |31| . The Fermi surface reconstruction re- 
sulting from the magnetic order on the UD side is not 
expected to lead to nodes, but may result in deep min- 
ima in the gap structure [321 133] . The increasing strength 
of the static order and low-frequency magnetic fluctua- 
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FIG. 3: Inhomogeneity in AAab(r) observed in an a; = 11% 
sample, a: Photograph of the x = 11% sample measured by 
SSS. The label of positions corresponds to the position where 
AAai,(r) was measured. b,c: l^\ab{T) by SSS at positions 
1-3 (b) and three different positions around position 3 each 
separated by 15/im (c). Dashed lines in (b) plot power-law 
fits with parameters given in Table |l] 



tions on the UD side [H [10] could enhance pair-breaking 
magnetic scattering |34j , giving rise to a power law depen- 
dence in AAafc(T) 131] that gets sharper with less doping 
on the UD side. On the OD side, deep gap minima may 
result from an anisotropic reduction of pairing strength 
as the doping moves further from the static magnetic 
order. Although pair-breaking may play some role and 
is one possible explanation for the spatial variation in 
the 11% sample, two facts suggest that our results are 
not dominated by pair-breaking processes from sample 
imperfection: the doping dependence we report is con- 
sistent with the low energy excitations measurements on 
annealed crystals [35 , and we observe flat AAab(r) at 
OptD. 

The third trend, the rapid increase Ps{T) of the UD 
and OptD when cooling through Tc, also agrees with the 
importance of magnetism. If the pairing is mediated by 
spin-fluctuations [4], forming superconductivity pushes 
the fluctuation spectrum to higher frequency, which fur- 
ther strengthens pairing, leading to a more rapid rise 
of Ps{T) than the standard BCS expression would give 
[34j . On the OD side, away from the magnetic order, the 
absence of low-frequency magnetic fluctuations may con- 
tribute to the slow rise of Ps{T) when cooling through T^. 
The slower rise at a; = 11% than the weak-coupling BCS 
result may hint that the two gaps have different Tc's. 

To conclude, by locally measuring \ab{T) and 
AAab(r), we observe systematic doping evolution of both 
the zero temperature value and the temperature depen- 
dence of Ps{T) in Ba(Fei_a;COa;)2As2 single crystals. Us- 
ing local scanning techniques, we reduce the error from 
sample inhomogeneity. The three systematic trends we 
observe on Ps{T) across the superconducting dome are 
consistent with the role of magnetism as a coexisting and 
competing order to the superconductivity as well as the 
pairing glue. 
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